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INTRODUCTION 


Cost model of phosphoric acid fuel cell powerplant includes two parts: a 

method for estimation of fuel cell system capital costs, and an economic 
analysis which determines the levelized annual cost of operating the system 
used in the capital cost program. 

Cost estimates are prepared for a given powerplant based on the equipment 
specifications discussed in the previous report of the performance model. 

Costs were estimated by determining the actual capacities of the equipment and 
the existing cost data. Current costs of these equipments in the form 
expected to be used were obtained from the references. Total module cost can 
be obtained by multiplying the equipment cost by the Direct Cost Factor (DCF), 
Indirect Cost Factor (ICF), and Contingency Factor (CF). 

The levelized annual cost of an investment is defined as the minimum 
constant net revenue required each year of the life of the project to cover 
all expenses, the cost of money, and the recovery of the initial investment. 
This is the capital investment analysis approach commonly used by electric 
uti lities. 

The cost model has been coded in Fortran programs with several input 
options. Mathematical formulation and program description will be discussed 
in this report. A sample problem will be presented to express the inputs and 
outputs. 
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I. SYSTEM DESCRIPTION 


As shown in Figure 1, methane which is circulated by compressor (C) is 
preheated by heat exchanger E-l prior to mixing it with the super heated steam 
which receives its heat by passing through heat exchanger E-9. Before entering 
the reformer, the methane steam mixture is heated via heat exchangers E-2 and 
E-3. Inside the reformer, methane is catalytically reformed by reaction with 
excess steam to produce carbon monoxide, carbon dioxide, and the desired pro- 
duct, hydrogen. The effluent from the reformer is cooled by flowing through 
heat exchanger E-2 before it enters the high temperature shift converter S-l. 
The function of the high temperature shift converter is to increase the hydro- 
gen concentration and to reduce the carbon monoxide concentration of the 
reformer gas effluent. The temperature of the effluent from the shift conver- 
ter S-l is then reduced by passing through heat excangers E-l, E-9 and E-6 
before entering the low temperature shift converter S-2. The low temperature 
shift converter further increases the hydrogen concentration by promoting the 
shift reaction at a lower operating temperature. The effluent from the low 
temperature shift converter then enters the fuel cell containing 1^, CO, 

CH^, CO^ and H^O. The fuel cell converts inputs of hydrogen and oxygen 
to DC power, water and heat. Oxygen is delivered to the fuel cell by air 
compressor A, which also provides air to the reformer burner. The spent fuel 
from the fuel cell anode goes to the burner after mixing with air supplied by 
compressor A. 

Before entering the burner, the mixture is preheated by the burner 
effluent via heat exchanger E-4. The spent fuel is then burned with whatever 
additional methane is needed to provide the thermal energy necessary for the 
reformer reaction. 
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Figure 1 Flow dl&graa of CSU designed PA1C systea 









Heat generated in the fuel cell is removed by heat exchangers E-7 and E-10. 
Heat from heat exchanger E-7 can then be utilized in industrial heat processing 
or space heating and cooling, while exchanger E-10 is used to preheat the water 
supplied by liquid separator Q to provide the necessary steam needed for the 
reforming process. The effluents from the burner and fuel cell cathode will 
have their water removed and separated by condenser. E-5 and liquid separator Q 
before allowing them to be exhausted to the atmosphere. 
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II. COST MATHEMATICAL MODEL 


2.1 Capital Investment 

Total module cost of a piece of equipment can be separated into two 
parts: F08 equipment cost and the working capital costs; the latter is 

related to the former. The relationship of total module cost and FOB 
equipment cost is shown in Figure 2, where the total module cost is obtained 
by multiplying the purchased equipment cost (FOB) by three factors: Direct 

Cost Factor (DCF), Indirect Cost Factor (ICF), and Contingency Factor (CF). 

The definitions of these are also shown in the figure. DCF and ICF of each 
equipment can be obtained from Refs. 3 and 4, where CF is the input option. 

The working capital cost is the difference of these two kinds of cost. 

All the costs were corrected by the Marshall and Swift cost index to be in 
constant mid-1981 dollars which is basic year used in the model. 

Equipment Cost 

There are several methods for estimating equipment cost. Three of them 
were used in the developed model for different components, which are power 
factor method, interpolation of true cost data, and unit-cost estimate. The 
fuel cell stack cost was estimated by unit-cost estimate method. For pumps 
and power inverter, linear interpolation was used to estimate the cost from 
tabulated data published by Exxon (Ref. 1). The power factor method was most 
used for the estimation of equipment cost in this model, which includes the 
reformer, the shift converters, the heat exchangers, the separator, and the 
compressors. 
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Figure 2 

GENERALIZED INVESTMENT COST ESTIMATING LOGIC (REF. 3) 


FOB Equipment 

Material Factor 
Labor Factor 

Direct M&L Cost 

Engineering 
Construction Overhead 

Bare Module Cost 

Contingency, etc. 

Fee 

Total Module Cost 


100 



XXX 

XX 


XXX 

— 


xx 

XX 


XXX 


XX 


X 


XXX 

<-* 


x.xx Direct Cost Factor (DCF) 
[Piping, Concrete, Steel, 
Instruments, Electrical, 
Insulation, Paint, & Labor] 


x.xx Indirect Cost Factor (ICF) 


x.xx Contingency Factor (CF) 


= 100 x DCF x ICF x CF 
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8riefly, the power factor method is 

C a ? 

S- =a i S + a 3 (1) 

where C = cost 

S = capacity 

a^, a^, and a^ are coefficients to be determined 

From (1) $n(| - a 3 ) = Qn a^ + a^ tin s (2) 

A linear regression on sample cost data will provide the values of a^, a^, and 
a^. Cost data have been obtained from the sources listed in the references. 

The linear interpolation algorithm is 

Y = YT(I-l) + [YT ( I )-YT( 1-1 )] [X-XT ( 1-1 ) ] / [XT ( I )-XT ( 1-1 ) ] (3) 

where Y is the cost of X capacity 

YT ( I ) is the listing cost of XT ( I ) listing capacity. 

The stack cost estimates were based on calculations of actual quantities of 
raw materials used to fabricate the components (unit-cost estimate). Current 
cost of raw materials, in the form expected to be used, were obtained from 
Chemical Marketing Report (Ref. 10) and Refs. 1 and 2. Fabrication costs were 
then determined by multiplying the material cost by a manufacturing cost 
factor, which was selected based on the production rate and the degree of 

* t 

automation envisioned for the manufacturing facility. The factor reflects 
manufacturing value added, including direct and supervisory labor plus other 
manufacturing burdens (e.g., maintenance and inventory costs). For example, 
the cost of catalyst (platinum) is 

CCP = (CPLxLCPxAAxNCELLxNS) x (1 + MCP) (4) 
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Energy Related (E): purchased power and fuel 

Non-Energy Related (NE): other variables and semi-variables 

Fixed Charges: depreciation, return-on-investment; income 

taxes, and local taxes and insurance. 

Those cost elements were first converted into a series of future cash flows 
(escalation allowed) which were then level ized to obtain a uniform annual cost 
series. This procedure is presented graphically in Figure 3. 

Level ized annual costs were determined from the following generalized 
relationship: 


N ( 1+i + e E) n 


LAC = IxFCR+E l XL CRF + NE [*= 

n =l (1+ Y ) n t " =1 

where FCR = fixed charge rate, and equal to 

CRF 


N (l + i + e NE) n 


and 


CRF 

(l + y) n Y 


(5) 


m, n 


B 


IW 


[1-t (DEP)-C] 


( 6 ) 


CRFm, ng: capital recovery factor for the after-tax cost of capital 
m and the economic life ng 


t: 

C: 

DEP: 


tax rate 

investment tax credit rate 

levelized depreciation factor (Sum of Years Digit) and 

z [ n , - 1/CRF ] 

L,I T m, n T J 


equal to 


n T (n T + 1) 


m 


(7) 


V 

m: 


tax depreciation life 

after tax cost of capital at the assumed inflation rate 
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+ 


I 


and 


total module cost in mid-1981 dollars, and equal to KmKeK (l + e. 




N*-No-L + 


Km 

L 

Ke 

K 

W 

e k 

N* 

No: 


cost-of-capital factor = e 0*418mL 
design and construction time 
escalation factor = e^'^^ e k + 
equipment cost 
working capital 

real capital cost escalation per year 
first year of commercial operation of the investment 
the year used as basis for the cost estimate k 
annual inflation rate 


E 

NE 

eE 

eNE 

Y 

n 

CRFr 


annual energy cost 

annual non-energy cost 

annual energy escalation 

annual non-energy escalation 

weighted cost of capital with inflation io 

project life 

capital recovery factor at y cost of capital and n years, which equal to 

(illilzi 

yU + y)" 


(8) 



CPL : 

cost of platinum, $/g 

LCP : 

2 

loading of platinum, 6/cm 

AA : 

2 

active area per cell, cm 

NCELL : 

number of cells per stack 

NS : 

number of stacks 

MCP : 

manufacturing factor for catalyst. 


The manufacturing cost factors used for estimating the cost of PAFC stack in 
this model were adopted from Ref. 1. More detailed description of this factor 
can be found in Ref. 4, pages 191-201. 


2.2 Levelized Annual Cost Analysis 

The levelized annual cost (LAC) of an investment is defined as the minimum 
constant net revenue required each year of the life of the project to cover all 
expenses, the cost of money, and the recovery of the initial investment. LAC 
is a comparative measure of both the fixed and variable costs associated with 
the investment, incurred at different times throughout the life of the project. 


The following formulations were taken principally from: NASA Documents 

dated April 1, 1979. Groundrules for Economic Analysis which also used in the 
study " Study of Component Technologies for Fuel Cell On-Site Integrated Energy 
Systems " , NASA CR— 165152 (December 1980), prepared by A. D. Little, Inc., for 
NASA Lewis Research Center. 

The computation of the levelized annual cost was accomplished by segre- 
gating annual costs into three categories, namely, energy related costs, 
non-energy related costs and fixed charges. The cost items grouped in each 
category were as follows: 
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Figure 3 


APPROACH TO LEVELIZED ANNUAL COST ANALYSIS 


PROJECT LIFE 



PROJECTED 
CASH FLOW 


I =Capital Investment 
NE = Non-energy Cost 
E = Energy Cost 




/ 


LEVELIZED 

ANNUAL 

COSTS 


K = Conversion Factor 
Defined as 


where: i = inflation 

e = real escalation 
n = year 

r = weighted cost of capital 


(1 + i+e) n 

ti <i+r » n 


FCR = Fixed Charge Rate with SYD 

depreciation and 10%- tax credit 
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III. COST COMPUTER MODEL 


3.1 Program 

There is one subroutine (RUN) in addition to the BLOCK DATA and MAIN 
programs in the cost computer model. The MAIN program estimates the capital 
investment of the PAFC powerplant, and calculates the levelized annual cost 
using the algorithm described in the previous chapter. The subroutine RLIN do 
the linear interpolation with two sets of input serial data and a specific 
capacity. The BLOCK DATA supplies the cost data tables, for the pump and the 
power inverter, from Ref. 1, and also the physical properties of the gases in 
the system. Table 1 shows the nomenclature of the variables. 

3.2 Program Operation 

The program input consists of a set of NAMELIST data which must be in a 
specified order. The first NAMELIST set is called INDEX and contains the 
Marshall and Swift cost index of the specified time. All the indices are 
obtained from Chemical Engineering magazine. 

The second set (CONST) has the constants used in the power factor method 
(Section 2.1). The general form used here is 

C = a L (S/a 2 ) a3 (9) 

where C is cost and S is capacity. The definitions of the constants for the 
equipment in this NAMELIST are listed in Table 2. 

The third set (FUCEC) contains the amount, the unit cost, and the manu- 
facturing cost factor of the material used in manufacturing the PAFC stack. 
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TABLE 1 


NOMENCLATURE OF COST COMPUTER MODEL 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 


Equipment Number and Unit for Estimating the Cost 


fuel cell stack 

reformer 

fuel compressor 

heat exchanger 

separator 

pump 

condenser 

high temperature shift converter 
low temperature shift converter 
power inverter 
air compressor 


kW 

MBtu/hr ejected 
brake HP 

transfer area ft^ 
g-mole water 
W 

gal /min water 
g-mole H2 
g-mole H2 
V 

f t^/min 


Cost of Fuel Cell Stack 


AA: 

NS: 

SV: 

CPL : 
CMRIN: 
NCELL : 
LCP: 
LESL : 
LEM: 
LBP 
CKW 
MCP : 
MESL : 
MEM: 
MBP: 
MCC 
MSH 
CCP: 
CGFP : 
CEM: 
CBP 
CCC 
CSH 
CGF 
CSC: 
CCPR: 
CMROT : 


index of raw material 


g/cm‘ 


active area per cell, cnr 
number of stacks 
operating voltage, V 
cost of platinum, $ 1 g 
Chemical Marketing Reporter 
number of cells per stack 
platinum loading, g/cm^ 
electrolyte support layers loading, 
electrolyte matrix loading, g/crcr 
bipolar plate loading, g/cm^ 
capacity of fuel cell stack, kW 
mfg. cost factor of catalyst 
mfg. cost factor of electrolyte support layers 
mfg. cost factor of electrolyte matrix 
mfg. cost factor of bipolar plate 
mfg. cost factor of cooling cartridge 
mfg. cost of factor of stack hardware 
cost of platinum (catalyst) 

cost of electrolyte support layers - graphite fiber paper 

cost of electrolyte matrix - silicon carbide fiber 

cost of bipolar plate - carbon/phenolic resin 

cost of cooling cartridge - carbon plate with copper tube grid 

cost of stock hardware - end plates, manifolding, tie rods 

unit cost of graphite fiber paper, £/g 

unit cost of silicon carbide fiber, $/ g 

unit cost of carbon/phenolic resin, $/g 

CMR index of data year 
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TABLE 1 (cont'd) 

NOMENCLATURE OF COST COMPUTER MODEL 


CC1: 

CC2: 

CPI: 

CP2: 

HCH4: 

HCO: 

HH2 : 

COST(I): 

CEQ(I.J): 

IN81: 

IN80: 

IN79: 

IN791: 

IN77 : 

IN75: 

IN68: 

IN67M: 

CH4: 

CO: 

H2: 

COMP : 

HE( J): 
SEPR: 

PUM: 

COND: 
HSHIF : 
LSHIF: 
AIRC: 


Cost of Other Equipments 

power conditioner voltage, V 
power conditioner cost, 2/kW 
pump power, W 
pump cost, % 

high heating value of methane, Cal/g-mole 

high heating value of carbon monoxide, Cal/g-mole 

high heating value of hydrogen, Cal/g-mole 

cost of equipment I, % 

capacity of equipment I number J 

Marshall and Swift index of mid-1981 

Marshall and Swift index of 1980 

Marshall and Swift index of 1979 

Marshall and Swift index of January 1979 

Marshall and Swift index of 1977 

Marshall and Swift index of 1975 

Marshall and Swift index of 1968 

Marshall and Swift index of mid-1967 

methane input, g-mole/hr 

carbon monoxide input, g-mole/hr 

hydrogen input, g-mole/hr 

brake hp of compressor, hp 

transfer area of heat exchanger number J, m^ 

amount of steam input in separator, g-mole/hr 

power of pump, hp 

inlet H 2 O flow rate of condenser 

inlet hydrogen flow rate of high temp, shift converter, g-mole/hr 
inlet hydrogen flow rate of low temp, shift converter, g-mole/hr 
inlet air flow rate, g-mole/hr 


Total Module Cost and Operation Cost 


DCF ( I ) : 
ICF(I): 
CF: 

CMAIN: 
CREPL: 
MTIME : 
WATER: 
CWAT : 
AVER: 
ENPU: 
AVHT : 
CENG: 


direct cost factor of equipment I 

indirect cost factor equipment I 

contingency factor of equipment 

maintenance cost of fuel cell system, S/kWh DC 

factor of capital cost for replacement 

times which replacement will occur for 20 years usage 

cooling water input, g-mole/hr 

cost of cooling water, $/ m^ 

mean factor of cooling water for recycle 

input fuel flow rate, g-mole/hr 

average heating value of input fuel, Btu/ft^ 

cost of energy fuel, 2/GJ 


< 
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TABLE 1 (cont'd) 

NOMENCLATURE OF COST COMPUTER MODEL 


CC: 

CD: 

CP: 

EK: 

ESC: 

FC: 

FD: 

FL: 

FP: 

L: 

NE: 

NSTAR: 

NT: 

NZERO: 

TAX: 

TAXL : 

TC: 

CAKE: 

CAKM: 

CAPIT: 

CEN: 

CN: 

CRFRE: 

CRFRK: 

CRFRT: 

DEP: 

FCL : 

FCR: 

R: 

RLAC: 

TLIN: 


Level i zed Annual Analysis 

cost of common equity 
cost of debt 

cost of preffered equity 

real capital cost escalation per year; i.e., rate of capital cost 

escalation, decimal 

ratio of common equity 

ratio of debt capital to total capital 

annual inflation rate 

ratio of preferred equity 

design and construction time, year 

economic life 

first full year of commercial operation of investment change above 
or below the rate of inflation 
tax depreciation life 
the year used as basic year 
tax rate 

state and local tax 
investment tax credit rate 
escalation factor 
cost-of-capital factor 
capital investment 
levelized energy cost 
non-energy cost 

capital recovery factor at R for economic life 

capital recovery factor at AK for energy in economic life 

capital recovery factor at R for tax depreciation life 

levelized depreciation factor for sum of years digits (SYD) 

levelized fixed charges 

fixed charge rate 

after tax cost of capital 

levelized annual cost 

levelized local tax and insurance 



The fourth set (INPUTS) consists of the input flow composition of fuel 
compressor, condenser, separator, high temperature and low temperature shift 
converters, the transfer area of each heat exchanger, and power needed in 
compressor and pump. 

The fifth set (FACTR) contains direct cost factor, indirect cost factor, 
and contingency factor of each equipment. 

The sixth and seventh sets (NENEG and ENG) include the amount and unit 
cost of fuel and utilities used in the system. The maintenance information is 
in NENEG. 

The last NAMELIST set (ECON) contains all the necessary data used for LAC 
analysis. 

All of the input variables are listed in Table 3, along with their units 
and numerical values in the sample run. 

3.3 Sample Problem 

The computer code described in the previous sections was used to estimate 

the equipment capital cost and the levelized annual cost of CSU designed PAFC 

powerplant (Figure 1). A 100 kW powerplant was considered here, which included 

2 

one fuel cell stack containing 200 cell plates with 1900 cm active area in 
each cell plate. The middle of year 1981 was chosen as the basic year for con- 
stant dollar estimation. 
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TABLE 2 


DEFINITIONS OF CONSTANTS 

IN NAMELIST 

CONST 



Constants 

Used in 

Equation 

Equipment 

a l 

a 2 

a 3 

Reformer 

Cl 

1 

C2 

Fuel Compressor 

C3 

1 

C4 

Heat Exchangers 

C5 

1 

C6 

Separator 

C7 

C8 

C9 

Pump 

CIO 

1 

Cll 

High Temperature Shift Converter 

C12 

C13 

C14 

Low Temperature Shift Converter 

C15 

C16 

C17 

Air Compressor 

C18 

1 

C19 
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TABLE 3 


— - - - - - - 


-INPUT- 

DATA OF SAMPLE PROBLEM - - 

NAMELIST 

Variable 

Sample 


Unit 

Definition 

Name 

Name 

Data 



. 

INDEX 

IN81 

696.9 



M.& S. index of mid 1981 

INDEX 

IN80 

659.6 



M.& S. index of 19 80 

INDEX 

IN791 

561 



M.& S. index of Jan. 1979 

INDEX 

IN77 

505.4 



M.& S. index of 1977 

INDEX 

IN75 

444 . 3 



M.& S. index of 1975 

INDEX 

IN6 8 

273 



M. & S. index of 1 968 

INDEX 

IN67M 

270 



M.& S. index of mid 1967 

CONST 

C 1 .... C 1 9 

7620 , 0 . 

85 

,514.55, 

0.82, 162. 106,0.6934, 1500, 



817200, 

0. 

64,104.4 

,0.5,900. ,4310,0.69,1320, 



4540,0. 

69 

,7,0.68 

constants listed in Table 2 

FUCEC 

AA 

1 900 


cm^ 

active area per cell 

FUCEC 

NS 

4 



number of stacks 

FUCEC 

SV 

133 


volt 

operating voltage in the 
stack 

FUCEC 

CPL 

16.75 


$/g 

cost of platinum( basic year) 

FUCEC 

CMRIN 

158.34 



CMR ( Chemical Marketing Report) 
index of raw material of basic 






year 

FUCEC 

NCELL 

200 



number of cells per stack 

FUCEC 

LCP 

0.00075 


g/cm 2 

loading of platinum 

FUCEC 

LESL 

0.024 


g/ cm 2 

loading of electrolyte support 
layers 

FUCEC 

LEM 

0.039 


g/cm 2 

loading of electrolyte matrix 

FUCEC 

LBP 

0.44 


g/cm 2 

loading of bipolar plate 

FUCEC 

CKW 

100 


. KW 

capacity of the fuel cell 

FUCEC 

MCP 

0.05 



mfg. cost factor of catalyst 

FUCEC 

MESL 

0.6 



mfg. cost factor of 
electrolyte support layers 

FUCEC 

MEM 

0 . 6 



mfg. cost factor of 
electrolyte matrix 

FUCEC 

MBP 

1 .5 



mfg. cost factor of bipolar 
plate 

FUCEC 

MCC 

1 .5 



mfg. cost factor of cooling 
plate 

FUCEC 

MSH 

1 .4 



mfg. cost factor of stack 
hardware 

FUCEC 

CGF 

0.066 


$/g 

unit cost of graphite fiber 
paper 

FUCEC 

CSC 

0.0176 


$/g 

unit cost of silicon carbide 
fiber 

FUCEC 

CCPR 

0.0009 


$/g 

unit cost of carbon/phenolic 
resin 

FUCEC 

CMROT 

198.66 



CMR index of data year 
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TABLE 3 

INPUT DATA OF SAMPLE PROBLEM 
continued 


NAMELIST Variable Sample Unit Definition 


Name 

Name 

Da ta 



INPUTS 

CH4 

172.6 

g-mole/ hr 

methane input flow rate 

INPUTS 

CO 

2.79 

g-mol e /hr 

carbon monoxide flow rate 

INPUTS 

H2 

867-63 

g-mole/hr 

hydrogen flow rate 

INPUTS 

COMP 

1.62 

hp 

brake hp of compressor 

INPUTS 

HE ( J ) 

0.3945, 

1.4024,1.5395,2.3735,1.4953,0.2,0.6418 
m^ transfer area of heat 

exchanger J 

INPUTS 

SEPR 

6820.63 

g-mole/hr 

input H20 flow rate of 
separator 

INPUTS 

PUM 

0.00226 

hp 

power of pump 

INPUTS 

COND 

132960 . 

37,51396 

g-mole/hr 

input H20 flow rate of 
condensers 

INPUTS 

HSHIF 

3708 . 6 

g-mole/hr 

input H2 flow rate of 
high temperature shift 
converter 

INPUTS 

LSHIF 

3925.62 

g-mole/hr 

input H2 flow rate of low 
temperature shift converter 

INPUTS 

AIRC 

24524 

g-mole/hr 

inlet air flow rate 

FACTR 

DCF ( I ) 

1.15,1. 
1 -75 

42,1.15,1 .35 

,1 .14,1 .75,1.16,1.15,1 .15,1.15 
direct cost factbr of 
equipment I 

FACTR 

ICF(I) 

1.14,1. 

1.14,1. 

28,1.14,1.1407,1.15,1.45,1.5086,1.14,1.14 
45 » indirect cost factor of 

equipment I 

FACTR 

CF 

0.2 


contingency factor of 
equipments 

NENEG 

CMAIN 

0 . 00065 

$/KW-h DC 

maintenance cost of system 

NENEG 

CREPL 

0.5 


factor of capital cost for 
replacement 

NENEG 

MTIME 

4 


times which replacement will 
occur for 20 yrs. 

NENEG 

WATER 

184356 

g-mole/hr 

cooling water flow rate 

NENEG 

CWAT 

.001316 

$/m3 

cost of cooling water 

NENEG 

AVER 

12 


mean factor of cooling water 
for recycle 

ENG 

ENPU 

1 405 . 1 6 

g-mole/hr 

input fuel flow rate 

ENG 

AVHT 

360242. 

6 Btu/ft3 

average heating value of input 
fuel 

ENG 

CENG 

6.29 

$/ G J 

cost of energy fuel 

ECON 

TAX 

0.48 


tax rate 

ECON 

TC 

0.1 


investment tax credit rate 

ECON 

ESC 

0.024 


escalation 

ECON 

CD 

0.03 


cost of debt 

ECON 

CP 

0.09 


cost of preferred equity 
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NAMELIST 

Name 


ECON 

ECON 

ECON 

ECON 

ECON 

ECON 

ECON 

ECON 

ECON 

ECON 

ECON 

ECON 


TABLE 3 

INPUT DATA OF SAMPLE PROBLEM 
continued 


Variable 

Name 

Sample 

Data 

Unit 

Definition 

CC 

0.09 


cost of common equity 

FD 

0.4 


ratio of debt capital to 
total capital 

FP 

0 


ratio of preferred equity 

FC 

0.6 


ratio of common equity 

TAXL 

0.02 


state and local tax 

FL 

0 


annual inflation rate 

NT 

20 


tax depreciation life 

NE 

20 


economic life 

L 

1 

year 

design and construction time 

EK 

0 


real capital cost escalation 
per year 

NSTAR 

1982 


first full year of commercial 
operation 

NZERO 

1981 


basic year 
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Figure 4 

SAMPLE INPUT DATA 


SINDEX IN81=696.9»IN80=659.6,IN791=561.,IN77=505.4»IN75=444.3,IN68=273. 
,IN67M=270. , 

SEND 

SCONST 01=7620. , C2= . 85 , C3 = 51 4 . 55 , 04 = . 82 , 05 = 162 . 1 06 , C6 = . 6 934 , 

07=150 0. ,08=81720 0 . ,C9=. 64, Cl 0 = 104.4 ,C11 =.5, 012 = 900 . 

013=4310. ,014= . 6 9, 015=1320. ,016=4540. ,017 =.6 9,018=7. ,C19=. 68, 

SEND 

SFUCEC AA=1900 . . NS=4 , SV=133 . , CPL=16 . 7 5 , CMRI N=1 58 . 34 , NCELL =200 , LCP= . 00 075 
,LESL=0.024,LEM=0.039,LBP=0.44,CKW=100. , MCP=0 . 05 , MESL =0 . 6 ,MEM=0 . 6 ,MBP=1 . 5 
,MCC=1.5,MSH=1.4,CGF=0.066,CSC=0.0176,CCPR=0.0009,CMROT=198.66, 

SEND 

S INPUTS CH4=172.6,C0=2.79,H2=867.63,C0MP=1 .62, 

HE=0. 3945, 1. 4024, 1. 5395, 2. 3735, 1.4953, 0.2,0. 6418,5 EPR=6820. 63, PUM=0. 00 226, 
COND=132960. 37,51396. ,HSHIF=3708.6,LSHIF=3925.62,AIRC=24524. , 

SEND 

SFACTR DCF=1. 15, 1.42, 1.15, 1.35, 1.14, 1.75, 1.16, 1.15, 1.15, 1.15, 1.75, 

I C F=l. 14, 1. 28, 1. 14, 1.407, 1.15, 1.45, 1.5086, 1.14, 1.14, 1.14,1. 45, CF=0. 2, 

SEND 

SNENEG CMAIN=0 .00065, OR EPL=0. 5, MTIME=4, WATER =184356., CWAT=Q .0013157, AVER =12. , 
SEND 

SENG ENPU=1405.16,AVHT=360242.64,CENG=6.29, 

SEND 

SECON TAX=0.48,TC=0.1,ESC=0.024,CD=0.03,CP=0.09,CC=0.09,FD=0.4,FP=0.,FC=0.6 
, T AXL = 0 . 02 , FL =0 . , NT=20 , NE=20 , L = 1 , EK = 0 . , NSTAR = 1 932 , NZERO = l 981 , 

SEND 
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Figure 5 

SAMPLE COMPUTER RUN 


SINDEX 

IN81 = 696.8999 

I N 8 0 = 659.5999 

IN791= 561.0 

IN77= 505.3999 

IN75 = 994.2998 

IN68 = 273.0 

IN67M- 270.0 

SEND 

SCONST 

Cl= 7620.0 

C2= 0.850 

C3= 514.5498 

C4= 0.820 

C5= 162.1060 

C6= 0.69340 

C7= 1500.0 

C8= 817200.0 

C9= 0.640 

C10= 104.40 

Cll= 0.50 

C12= 900.0 

C13= 4310.0 

C 1 4 - 0.690 

C15= 1320.0 

Cl 6 = 4540 . 0 

C17= 0.690 

C18= 7.0 

Cl 9= 0.6799999 

SEND 

SFUCEC 

A A - 1900.0 

NS = 4 

SV= 133.0 

CPL = 16.750 

CMRIN- 158.340 

NCEl L = 200 

LCP= 0 . 7499999E-03 

LESL= 0.240E-01 

IEM = 0.390E-01 

LBP= 0.440 

CKW= 100.0 

MCP = 0.50E-01 

MESL = 0.60 

MEM= 0.60 

M3P = 1.50 

MCC = 1.50 

r,SH= 1.4 0 

CGF = 0.6599998E-01 

CSC= 0. 1760 E- 01 

CCPR = 0 .8999999E-03 

CMROT= 198.660 

SEND 

SINPUTS 

CH4= 172.60 

CO= 2.790 

H2 - 867.6299 

COMP = 1.620 

HE= 0.39450, 1.402399, 1.539499, 2.37350, 1.495299, 0.20, 
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Figure 5 (cont'd) 
SAMPLE COMPUTER RUN 


SEPR= 6820.629 

PUM = 0.2260E-02 

C0ND = 132960.3, 51396.0 

HSH I F= 3708.60 

LSHIF= 3925.620 

AIRC= 29529.0 

& END 

8FACTR 

ICF= 1.139999, 1.280, 1.139999, 1.9070, 1.150, 1.950, 1.508599, 3X1.139999 
1.950 

DCF= 1.150, 1.919999, 1.150, 1.399999, 1.139999, 1.750, 1.160, 3*1.150, 1.750 
CF= 0.20 
& END 
SNENEG 

CMAIN= 0 . 6999998E-03 
CREPL = 0.50 
MTIME= 9 
WATER= 189356.0 
CWAT = 0 . 131570E-02 
AVER= 12.0 
XEND 
& ENG 

ENPU = 1905.160 
AVHT= 360292.6 
CENG= 6.290 
& END 

COST ANALYSIS FOR 100KW FUEL CELL SYSTEM 

MID-1981 MONEY 
100% LOAD FACTOR 


EQUIPMENT 

CAPITAL COST ( F. 0 . B ) 

PERCENTAGE 

COST ( 

11 = 

0.28001E 

05 

99.80 

COST ( 

2) = 

0 .85823E 

09 

13.73 

COST ( 

3) = 

0 . 19509E 

09 

3.12 

COST ( 

9) = 

0 .76818E 

09 

12.29 

COST ( 

51 = 

0 . 96691E 

02 

0 . 15 

COST ( 

6 ) = 

0.52895E 

03 

0.85 

COST ( 

7) = 

0 . 19186E 

09 

2.27 

COST ( 

8) = 

0.11188E 

09 

1.79 

COST ( 

9 ) = 

0 . 16969E 

09 

2.63 

COSTC 1 0 ) = 

0 . 10533E 

05 

16.85 

COSTC 11 ) = 

0.93990E 

03 

1.50 


TOTAL CAPITAL COST(F.O.B) 0.62997E 05 
TOTAL WORKING CAPITAL COST 0.36873E 05 
ANNUAL OSM 0.83828E 09 

ANNUAL ENERGY COST INYEAR J=0 0.61990E 05 


SECON 

TAX= 0.980 
TC= 0.9999996E-01 
ESC= 0.290E-01 
CD= 0.30E-01 
CP= 0 .8999997E-01 
CC= 0.8999997E-01 
FD= 0.90 
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Figure 5 (cont'd) 
SAMPLE COMPUTER RUN 


FP= 0.0 

FC= 0.60 

TAXL = 0.20E-01 

FL= 0.0 

NT = 20 

NE= 20 

L= 1 

EK= 0.0 

N5TAR= 1982 

NZERO= 1981 

SEND 

INFORMATION OF ECONOMIC FACTOR: 

LEVELIZED DEPRECIATION FACTOR (SYD) 0.67699 
FIXED CHARGE RATE 0.09791 

CAPITAL RECOVERY FACTOR OF ECOMONIC LIFE 0.08718 

CAPITAL RECOVERY FACTOR OF TAX DEPRECIATION LIFE 0.08718 


LEVELIZED FIXED CHARGES 0.98846E OA 
LEVELIZED ENERGY COST Q.7608AE 05 

TOTAL LEVELIZED COST 0.97380E 05 
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The following are the summary of the results: 


1. Equipment Capital Cost (FOB) - in mid-1981 money 


Equipment 

Cost (FOB)-J 

Percentage of 
Total FOB 

fuel cell module 

28001 

44.8 

reformer 

8582 

13.7 

fuel compressor 

1951 

3.1 

heat exchangers 

7682 

12.3 

separator 

97 

0.2 

pump 

528 

0.9 

condenser 

1419 

2.3 

high temperature shift converter 

1119 

1.8 

low temperature shift converter 

1646 

2.6 

power inverter 

10535 

16.8 

air compressor 

939 

1.5 

total 

62497 

100.0 


2. Total Working Cost 

Total Working Cost = total module cost - total FOB cost (Figure 2) 
36873 = 99370 - 62497 


3. Levelized Annual Analysis 


annual operation and maintenance 8383 
levelized local tax and insurance 3028 
levelized energy cost 76084 
levelized fixed charges 9885 


total levelized annual cost 97380 


The required CPU time to run this sample problem is less than 0.01 minute on 
IBM/370. 
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LISTING OF THE COST COMPUTER MODEL 
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ro 

00 


0000100 

0000200 

0000300 

0000400 

0000500 

0000600 

0000700 

0000800 

0000900 

0001000 

0001100 

0001200 

0001500 

0001600 

0001700 

0001920 

0001940 

0002120 

0002140 

0002200 

0002400 

0002500 

0002600 

0002700 

0002800 

0002900 

0003000 

0003100 

0003200 

0003300 

0003400 

0003500 

0003600 

0003700 

0003800 

0003900 

0004000 

0004100 

0004200 

0004300 

0004400 

0004500 

0004600 

0004700 

0004800 

0004900 

0005000 

0005100 

0005200 

0005300 

0005400 

0005500 

0005600 

0005700 


CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C THIS PROGRAM IS TO CALCULATE GENERALIZED INVESTMENT COST ESTIMATING LOGIC x 
C WHICH IS RECOMMENDED BY K. M. GUTHRIE, "PROCESS PLANT ESTIMATING, EVALUATIONx 
C , AND CONTROL" x 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
BLOCK DATA 

REAL CC1C20),CC2(20),CP1(20),CP2(20) 

COMMON/DAT A/CC 1 ,CC2, CPI, CP2 , HCH4 , HCO , HH2 
C CCl: POWER CONDITION VOLT (VOLT) 

C CC2 : POWER CONDITION COST ($/KW) 

C CP1-- PUMP POWER (WATT) 

C CP2 : PUMP COST ( $ ) 

C HCH4 : HIGH HEAT VALUE OF CH4 (CAL/G-MOLE) 

C HCO: HIGH HEATING VALUE OF CO (CAL/G-MOLE) 

C HH2 : HIGH HEATING VALUE OF H2 (CAL/G-MOLE) 

DATA CC1/50. ,164., 203., 248., 304., 36 6., 433., 528., 657 
DATA CC1(12)/2810./,CC1( 13)/100000 0000 ./ 

DATA CC2/200 .,160. ,150. ,140. ,130. ,120. ,110. ,100. ,90 
DATA CC2 ( 1 3 )/50 . / 

DATA CP 1/0 . ,6150 0 . ,264000 . , 6 150 0 0 . / , CP2/50 0 . , 6 7 0 0 . , 

DATA HCH4/2 12800 ./,HCO/67636 ./.HH2/68317 ./ 

END 

REAL ICF(11),DCF(11),CEQ(20,10),COST(20),CC1(20),CC2(20),CP1(20) - 

1 , CP2 ( 20 ) , LCP,LESL,LEM,LBP,MCP,MESL,MEM,MBP,MCC,M5H - 

2, LSHIF, IN81,IN80,IN7 91,IN68, IN6 7M, IN7 7,IN7 5 
DIMENSION P(11),HE(7),C0ND(2) 

COMMON/DATA/ CCl , CC2 , CP 1 , CP2 , HCH4 , HCO , HH2 

NAMEL1ST/FUCEC/ A A , NS , SV , CP L , CMRIN , NCELL , LCP , LESL , LEM, LBP , CKW 
1,MCP,MESL,MEM,MBP,MCC,MSH,CGF,CSC,CCPR,CMR0T 
NAMELIST/ INPUTS/ CH4 , CO , H2 , COMP , HE , SEPR , PUM , COND , HSHI F , LSHIF , AIRC 
NAMELIST/INDEX/ IN81 , 1 N80 , IN7 9 1 , IN77 , IN75 , IN6 8 , IN6 7M 
NAMELIST/FACTR/ ICF,DCF,CF 

NAMELIST/NENEG/ CMAI N , CREPL , MT IME, WATER , CWAT , AVER 
NAMELIST/ENG/ ENPU , AVHT , CENG 

NAMELIST/ ECON/ T AX , TC , ESC , CD , CP , CC , FD , FP , FC , TAXL , FL , NT , NE, L , EK , 

INSTAR, NZERO 

NAMELIST/CONST/ Cl , C2 , C3 , C4 , C5 , C6 , C7 , C8 , C 9 , Cl 0 , Cl 1 , Cl 2 , Cl 3 , C14 , 
1C15,C16,C17,C18,C19 
C 

CXXXXXXXXXXXXXXkXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C EQUIPMENT NO. AND UNIT FOR CALCULATING COST 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C is FUEL CELL, KW 
C 2- REFORMER, MBTU/HR EJECTED 
C 3: COMPRESSOR(GAS) , BRAKE HP 
C 4: HEAT EXCHANGER, TRANSFER AREA FTXX2 
C 5: SEPARATOR, G-MOLE H20(L) 

C 6 : PUMP, WATTS 
C 7: CONDENSER, GAL. /MIN. H20(L) 

C 8 : SHIFT CONVERTER(HIGH TEMPERATURE), MOLES H2 
C 9: SHIFT CONVERTER( LOW TEMPERATURE), MOLES H2 
C 10: POWER INVERTER, SYSTEM VOLT 
C 11= AIR COMPRESSOR (BLOWER), FTXX3/MIN. 

C 


. , 920. ,1560./ 
.,80,70. ,60./ 
32000 ., 95400 ./ 



no 


0005800 

0005900 

0006000 

0006100 

0006200 

0006300 

0006900 

0006500 

0006600 

0006700 

0006800 

0006900 

0007000 

0007100 

0007200 

0007300 

0007900 

0007500 

0007600 

0007700 

0007800 

0007900 

0008000 

0008100 

0008200 

0008300 

0008900 

0008500 

0008600 

0008700 

0008800 

0008900 

0009000 

0009100 

0009200 

0009300 

0009900 

0009500 

0009600 

0009700 

0009800 

0009900 

0010000 

0010100 

0010200 

0010300 

0010900 

0010500 

0010600 

0010700 

0010800 

0010900 

0011000 

0011100 


CXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C DEFINITION: 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C COST ( I ) : COST OF EQUIPMENT I 

C CEQ( I , J ) :■ CAPACITY OF EQUIPMENT I NO.J (ACCORDING TO THE COST ESTIMAT 

C 

C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C INPUT FUNCTIONS FOR CALCULATING COST OF EACH EQUIPMENT 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C 

C BASIS :MID~1981 MONEY 
C 100X LOAD FACTOR 

C 

F2(S)=C1X(S)*XC2XIN81/IN68 

F3(S)=C3X(S)XXC9XIN81/IN68 

F9(S)=C5XSXXC6XIN81/IN791 

F5(S)=C7X(S/C8)XXC9XIN81/IN77 

F7(S)=C10X(S)XXC11XIN81/IN67M 

F8(S)=C12X(S/C13)XXC19XIN81/IN77 

F9(S)=C15X(S/C16)XXC17XIN81/IN77 

F11(S)=C18XSXXC19XIN81/IN68 

C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C READ IN THE MARSHALL AND SWIFT INDEX 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C IN81 : INDEX OF MID 1-981 
C IN80: INDEX OF 1980 
C IN7 9 ! 4HDEX OF 1979 
C IN791 : INDEX OF 1979 JAN. 

C IN7 7 : INDEX OF 1977 
C IN75 : INDEX OF 1975 
C IN68 : INDEX OF 1968 
C IN67M: INDEX OF MID. 1967 
C 

READC5, INDEX) 

WRITE(6, INDEX) 

READC5, CONST) 

WRITEC6, CONST) 

C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C CAL. THE COST OF FUEL CELL 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C 

C INPUT: 

C AA: ACTIVE AREA PER CELL (CMXX2) 

C NS= NUMBER OF STACKS 
C SV: STACK VOLTAGE( VOLT ) 

C CPL : COST OF PI ATINUM( $/G) — BASED ON BASIC YEAR 

C CMRIN : CMR(CHEMICAL MARKETING REPORTER) INDEX OF RAW MATERITAL OF BASI 
C NCELL : NUMBER OF CELLS PER STACK 
C LCP : LOADING OF PIATINUM(G/CMXX2) 

C L ESL : LOADING OF ELECTROLYTE SUPPORT LAYERS ( G/CMXX2 ) 

C LEM: LOADING OF ELECTROLYTE MATRIXC G/CMXX2 ) 



oo 


o 


0011200 C LBP: LOADING OF BIPOLAR PL ATE( G/CM**2 ) 

0011300 C CKW • CAPACITY OF THE FUEL CELL (KW) 

0011400 C MCP= MFG . COST FACTOR OF CATALYST 

0011500 C MESL:MFG. COST FACTOR OF ELECTROLYTE SUPPORT LAYERS 

0011600 C MEM: MFG. COST FACTOR OF ELECTROLYTE MATRIX 

0011700 C MBP: MFG. COST FACTOR OF BIPOLAR PLATE 

0011800 C MCC: MFG. COST FACTOR OF COOLING CARTRIDGE 

0011900 C MSH : MFG. COST FACTOR OF STACK HARDWARE 

0012000 C CCP : COST OF CATALYST-- PLATTINUM 

0012100 C CGFP : COST OF ELECTRODE SUPPORT LAYERS-- GRAPHITE FIBER PAPER 
0012200 C CEM= COST OF ELECTROLYTE MATRIX-- SILICON CARBIDE FIBER 
0012300 C CBP : COST OF BIPOLAR PLATE-- CARBON/PHENOLIC RESIN 

0012400 C CCC : COST OF COOLING CARTRIDGE-- CARBON PLATE WITH COPPER TUBE GRID 

0012500 C CSH : COST OF STACK HARDWARE-- END PLATES, MANIFOLDING, TIE RODS 

0012600 C CGF : UNIT COST OF GRAPHITE FIBER PAPER, $/G 

0012700 C CSC: UNIT COST OF SILICON CARBIDE FIBER, $/G 

0012800 C CCPR : UNIT COST OF CARBON/PHENOLIC RESIN, $/G 

0012900 C CMROT: CMR INDEX OF DATA YEAR 

0013000 C 

0013100 READ( 5 , FUCEC ) 

0013200 WRITEC6, FUCEC) 

0013300 CCP=(CPL*LCP*AA*NCELL#NS)#(1.+MCP) 

0013400 CGFP=(CGF*LESL*AAXNCELL*NSXCMRIN/CMR0T)X(1.+MESL) 

0013500 CEM=(CSCXLEMXAAXNCELLXNS*CMRIN/CMR0T)X(1 . +MEM) 

0013600 CBP=CCCPRXLBPXAAXNCELLXNSXCMRIN/CMR0T)X( 1 ,+MBP) 

0013700 C ASSUME THE RAW MATERITAL COST OF COOLING CARTRIDGE AND STACK HARDWARE 
0013800 C IS THE SAME AS BIPOLAR PLATE 
0013900 CCC= CBP/( 1 . +MBP ) X ( 1 . +MCC) 

0014000 CSH=CBP/(l.+MBP)x(l.+MSH) 

0014100 COST(l)=CCP+CGFP+CEM+CBP+CCC+CSH 

0014200 C 

0014300 CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

0014400 C INPUT THE CAPACITY OF EACH EQUIPMENT AND CALCULATE THE COST 

0014500 cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

0014600 C CH4: CH4 INPUT , G-MOL E/HR 

0014700 C CO: CO INPUT , G-MOL E/HR 

0014800 C H2 : H2 INPUT, G-MOLE/HR 

0014900 C COMP: BRAKE HP OF COMPRESSOR 

0015000 C HE= TRANSFER AREA OF HEAT EXCHANGER, Mxx2 

0015100 C SEPR: AMOUE r OF H20 INTO SEPARATOR , G-MOL E/HR 

0015200 C PUM: POWER OF PUMP, HP 

0015300 C COND: AMOUNT OF H20 INTO CONDENSER , G-MOLE/HR 

0015400 C HSHIF: AMOUNT OF H2 INTO HIGH TEMP. SHIFT CONVERTER, G-MOLE/HR 
0015500 C L SHI F : AMOUNT OF H2 INTO LOW TEMP. SHIFT CONVERTER, G-MOLE/HR 
0015600 C AIRC : INLET AIR, G-MOLE/HR 
0015700 C 

0015800 READC5, INPUTS) 

0015900 WRITE(6, INPUTS) 

0016000 CEQC2,l)=(CH4XHCH4+C0XHC0+H2XHH2)X3.97E-3/l.E+6 

0016100 C0ST(2)=F2CCEQ(2,1)) 

0016200 CEQ(3,1)=C0MP 

0016300 COST(3)=F3(CEQ(3,l)) 

0016400 COST ( 4 ) =0 . 


0016500 DO 1 K=l,7 




0016600 

0016700 

0016800 

0016900 

0017000 

0017100 

0017200 

0017300 

0017600 

0017500 

0017600 

0017700 

0017800 

0017900 

0018000 

0018100 

0018200 

0018300 

0018600 

0018500 

0018600 

0018700 

0018800 

0018900 

0019000 

0019100 

0019200 

0019300 

0019600 

0019500 

0019600 

0019700 

0019800 

0019900 

0020000 

0020100 

0020200 

0020300 

0020600 

0020500 

0020600 

0020700 

0020800 

0020900 

0021000 

0021100 

0021200 

0021300 

0021600 

0021500 

0021600 

0021700 

0021800 

0021900 


‘ CEQ (6, K)-HE(K)/. 3068**2 

1 C0ST(6)=C0ST(6)+F6(CEQ(6,K)> 

CEQ( 5,1) = SEPR 
C0ST(5)=F5(C£Q(5,1')) 

CEQ(6, 1)=PUM*765.7 

CALL RLIN(6,CP1,CP2,CEQ(6,1),C0ST(6)) 

COST(6)=COST(6)*IN81/IN80 
CEQ(7,l)=COND(l)*18./1000./3. 785/60. 

CEQ(7,2)=COND(2)*18./1000 ./3. 785/60 . 
COST(7)=F7(CEQ(7,l))+F7(CEQ(7,2)) 

CEQ ( 8 » 1 ) =HSHI F 
C0ST(8)=F8(CEQ(8,1)) 

CEQ(9,1)=LSHIF 

COST(9)=F9(CEQ(9,l)) 

CEQ( 10,1 ) =SV*NS 

CALL RLIN(13,CC1,CC2, CEQ (10,1) .COST (10)) 

C05T(10)=(COST(10)*IN81/IN80)*CKW 

CEQ( 11 ,1 )=AIRC/653.6*10 .73*298 .*1 .8/16.7/1 . 06/6 0 . 

C0ST(11)=F11(CEQ(11»1)) 

CAK=0. 

DO 2 K=1 , 1 1 

2 CAK=CAK+COST(K) 

DO 3 K=l,ll 

3 P(K)=COST(K)/CAK*100. 

C*************************************************************** ******** 
C INPUT DIRECT AND INDIRECT COST FACTORS 

C* ********************************************************************** 
C DCF( I ) •" DIRECTOR COST FACTOR OF EQUIPMENT I 
C ICF ( I ) • INDIRECT COST FACTOR OF EQUIPMENT I 
C CF : CONTINGENCY FACTOR OF EQUIPMENT 
READ( 5 , FACTR ) 

WRITE( 6 , FACTR ) 

DO 6 K=1 , 11 

6 CAW=CAW+COST(K)*(DCF(K)*ICF(K)-l. ) 

CAW=CAW*(CF+1 . ) 

C************************************** (******************************** 
C INPUT THE OPERATING AND MAINTENANCE COSTS (NONENERGY) 

C****** ********************************************************** ******* 
C CMAIN J MAINTENANCE COST OF FUEL CELL, $/KWH DC OUTPUT 
C CREPL : FACTOR OF CAPITAL COST FOR REPLANCEMENT 
C MTIME: TIMES WHICH REPLACEMENT WILL OCCUR FOR 20 YRS USAGE 
C WATER: INPUT COOLING WATER, G-MOLE/HR 
C CWAT COOLING WATER COST, $/M**3 
C AVER: MEAN FACTOR OF COOLING WATER FOR RECYCLE 
C 

READ (5.NENEG) 

WRITE ( 6 , NENEG ) 

OANDM=CKW*CMAIN*26.*365.+CAK*CREPL/MTIME+WATER 
1*18-/1000000. *CWAT*26 . *AVER 

C*********************************************************************** 
C INPUT THE ENERGY COST THEN CAL. ENERGY OPERATING COST 

C*********************************************************************** 
C ENPU : TOTAL INPUT FUEL, G-MOLE/HR 

C AVHT : AVERAGE HEATING VALUE OF INPUT FUEL, BTU/FT**3 



0022000 C CENG* COST OF ENERGY FUEL, $/GJ 
0022100 C 

0022200 READ( 5 , ENG) 

0022300 WRITE(6 , ENG) 

0022400 P0=ENPU/453.6*AVHT/1000000.*CENG#24.x365. 

0022500 C WRITE THE RESULTS 
0022600 WRITEC6 ,103) 

0022700 WRITE(6,ldl) ( ( KK , COST ( KK ) , P ( KK ) ) , KK=1 , 1 1 ) 

0022800 WRITE(6,102) CAK , CAW, OANDM , PO 

0022900 C ... 

0023000C#####8####8#######ft####«######ft########8#8#ft###################8#8###ft# 
0023100 C PERFORM THE ECONOMIC CALCULATION AND A CASH FLOW ANALYSIS 
0 02320 0 Cftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftiftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftft 
0023300 C 

0023400 CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
0023500 C INPUT THE ECONOMIC ANALYSIS FACTOR 

0023600 cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
0023700 C TAX s TAX RATE 

0023800 C TC: INVESTMENT TAX CREDIT RATE 
0023900 C ESC: ESCALATION, DECIMAL 
0024000 C CD: COST OF DEBT 
0024100 C CP: COST OF PREFERRED EQUITY 
0024200 C CC: COST OF COMMON EQUITY 

0024300 C FD: RATIO OF DEBT CAPITAL TO TOTAL CAPITAL 

0024400 C FP= RATIO OF PREFERRED EQUITY 

0024500 C FC: RATIO OF COMMON EQUITY 

0024600 C T AXL : STATE AND LOCAL TAX 

0024700 C FL= ANNUAL INFLATION RATE 

0024800 C NT: TAX DEPRECITION LIFE 

0024900 C NE: ECOMONIC LIFE 

0025000 C L= DESIGN AND CONSTRUCTION TIME, IN YEAR 

0025100 C EK *■ REAL CAPITAL COST ESCALATION PER YEAR, I. E., THE RATE OF CAPITAL 
0025200 C CHANGE ABOVE OR BELOW THE RATE OF INFLATION 

0025300 C NSTAR : FIRST FULL YEAR OF COMMERCIAL OPERATION OF THE INVESTMENT 
0025400 C NZERO: THE YEAR USED AS BASIS FOR THE COST ESTIMATE 
0025500 C 

0025600 READ(5, ECON) 

0025700 WRITE! 6 , ECON ) 

0025800 C R: AFTER TAX COST OF CAPITAL 

0025900 R=( 1 . rAX+TAXL))*FDXCD+FPXCP+FC*CC+FL.x(l.-(TAX+TAXL)*FD) 

0026000 C CAKM : COST-OF-CAPITAL FACTOR 
0026100 CAKM = EXP( ,418*R*L ) 

0026200 C CAKE: ESCALATION FACTOR 
0026300 CAKE=EXP( .562*(EK+FL)*L) 

0026400 C CAPIT: CAPITAL INVESTMENT 

0026500 CAPIT=CAKMXCAKE*CAKX(1 . +EK+FL )**(NSTAR-NZERO-L l+CAW 

0026600 C TLIN : LEVELIZED LOCAL TAX AND INSURANCE 

0026700 TLIN=0 . 03XCAPIT 

0026800 C CN*- NON-ENERGY COST 

0026900 CN=OANDM+TL IN 

0027000 C CRFRE : CAPITAL RECOVERY FACTOR AT R FOR ECOMONIC LIFE 
0027100 C3=l. 

0027200 C4=0. 

0027300 DO 5 1=1, NE 
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- C3=C3/(1.+R) 

C4=C4+C3 

5 CONTINUE 
CRFRE=1 . /C4 

C CRFRT : CAPITAL RECOVERY FACTOR AT R FOR TAX DEPRECIATION LIFE 
Dl = l. 

D2 = 0 . 

DO 6 1=1, NT 
D1=D1/(1 .+R) 

D2=D2+D1 

6 CONTINUE 
CRFRT = 1 . /D2 

C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C CALCULATION ANNUAL COST OF ENERGY ( VARY AT A CONST. ANNUAL RATE 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C CRFRK : CAPITAL RECOVERY FACTOR AT AK FOR ENERGY IN ECOMONIC LIFE 
AK= (l.+R)/(l.+ESC+FL)-l. 

G1 = 1 . 

G2 = 0. 

DO 7 J = 1 » NE 
G1=G1/(1 .+AK) 

G2=G2+G1 

7 CONTINUE 
CRFRK=1 ./G2 

C CEN: LEVELIZED ENERGY COST 
CEN=P0XCRFRE/CRFRK 

C DEP= LEVELIZED DEPRECIATION FACTOR FOR SUM OF YEARS DIGITS (SYD) 

DEP=2. xc NT- 1 ./CRFRT >/<NTX<NT+l . )XR) 

C FCR: FIXED CHARGE RATE 

FCR=(CRFRE/(1 .-(TAX+TAXL ) ) ) x ( 1 . -( TAX+TAXL ) XDEP-TC) 

C RL AC : LEVELIZED ANNUAL COST 
RLAC=CAPITXFCR+CN+CEN 
C FCL : LEVELIZED FIXED CHARGES 
FCL=CAPITXFCR 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C WRITE THE RESULTS 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
WRITEC6 , 104 ) 

WRITER, 106) DEP, FCR, CRFRE, CRFRT 
WRITE(6,105) FCL , CEN, RL AC 
C 

101 FORMAT (IX, ’COST( » ,12, * ) = ' , E13 . 5 , 10X, F5 . 2 ) 

102 F0RMAT(//1X, 'TOTAL CAPITAL COST ( F . 0 . B ) ' , El 3 . 5/1X , ' TOT AL WORKING CA- 
1PITAL COST’ .E13.5/1X, 'ANNUAL OJM ’ , El 3 . 5/ IX , ’ ANNUAL ENERGY COST IN- 
2YEAR J=0’ ,E13.5//) 

103 FORMAT (/ ' COST ANALYSIS FOR 100KW FUEL CELL SYSTEM'//1X, 'MID-198- 
11 MONEY '/IX, ’ 100/C LOAD FACTOR '//IX, ' EQUIPMENT CAPITAL COST(F.O.B) - 
2 PERCENTAGE') 

104 FORMATC/1X, ' INFORMATION OF ECONOMIC FACTOR:'/) 

105 FORMAT ( ' LEVELIZED FIXED CHARGES ’,E13.5//' LEVELIZED ENERGY COST - 
1 ' , El 3 . 5/// ' TOTAL LEVELIZED COST ’,E13.5) 

106 FORMAT < IX , ' LEVELIZED DEPRECIATION FACTOR (SYD) ',F10.5/ 

1 IX,’ FIXED CHARGE RATE ',F10.5/ 



3032800 

3032900 

3033000 

3033100 

3033200 

3033300 

3033900 

3033500 

3033600 

0033700 

3033800 

0033900 

0039000 

0039100 

0039200 

0039300 

0039900 

0039500 

0039600 

0039700 


2 IX,’ CAPITAL RECOVERY FACTOR OF ECOMONIC LIFE ’,F10.5/ 

3 IX,’ CAPITAL RECOVERY FACTOR OF TAX DEPRECIATION LIFE ’ 

9 , F10 . 5// ) 

C 

STOP 

END 

SUBROUTINE RL IN ( N , XT , YT, X, ANS ) 

C THIS SUBROUTINE IS TO CAL. LINEAR INTERPOLATION. 

C THE ALGORITHM REQUIRES XT VECTOR TO BE IN ASCENDING ORDER 

DIMENSION XT (20), YT (20) 

1=2 

IF(X. LE.XT(l)) GO TO 20 
I =N 

IF(X.GE.XT(N)) GO TO 20 
DO 10 1=2, N 

IF(X.LE.XTd)) GO TO 20 
10 CONTINUE 

20 ANS=YT(I-1)+(YT(I)-YT(I-1))/(XT(I)-XT(I-1))*CX-XT(I-1>) 

RETURN 

END 
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